HIGH SPEED EXCITED MULTI-SOLITONS IN NONLINEAR SCHRODINGER 

EQUATIONS 

RAPHAEL COTE AND STEFAN LE COZ 

>—^ ' Abstract. We consider the nonlinear Schrodinger equation in R'' 

^^ ■ idtu + Am + /(n) = 0. 






For d ^ 2, this equation admits travcUing wave solutions of the form e"^*$(x) (up to a Galilean 
transformation), where # is a fixed profile, solution to — A<J> + ui^ = /($), but not the ground 
state. This kind of profiles are called excited states. In this paper, we construct solutions to 
^+ ' NLS behaving like a sum of N excited states which spread up quickly as time grows (which we 

^Nl , call multi-solitons). We also show that if the flow around one of these excited states is linearly 

unstable, then the multi-soliton is not unique, and is unstable. 



1. Introduction 



< 

Cd , Setting of the problem. We consider tlie nonlinear Schrodinger equation 

(NLS) iut + Au + f{u) = 

where w : K x R'* ^ C and / : C ^ C is defined for any z e C hy f{z) = 
CN ; ggC°([0,+oo),R)nCH(0,+oo),R). 



Equation (NLS) admits special travelling wave solutions called solitons: given a frequency 
Wo > 0, an initial phase 70 € K, initial position and speed xq, uq G K'' and a solution $0 G i?^(K'') 



O ■ of 



(1) ^A$o+^o$o-/($o)=0, 

a soliton solution of ( NLS| ) travelling on the line a; = xq + v^t is given by 



cn 

o 

5 ■ (2) i?*o,-o,7o,«o,xo(t,a;):=<i>o(x-i;ot-xo)e'(^''''-"-^l"°l'*+"°*+^°). 

Among solutions of (|l|), it is common to distinguish between ground states, and excited states. A 
ground state (or least energy solution) minimizes among all solutions of (|l|) the action Sq, defined 
^ ■ for we H^{W^) by 

Cd : So{v) := -\\yv\\l2(jf.a) + Yl|w||i2(Rd) - j ^ F{v)dx, 

where F{z) := J^.'^ g{s^)sds for all z g C An excited state is a solution to (Q) which is not a 
ground state. In general, we shall refer to any solution of (|^) as bound state. We also mention the 
existence of a particular type of excited states, the vortices. A vortex is a special solution of (|l|) 
which is non-trivially complex-valued, i.e. with a non-zero angular momentum. Vortices can be 
constructed following the ansatz described by Lions in [Q . We shall sometimes abuse terminology 
and call ground state (resp. excited state) a soliton build with a ground state (resp. an excited 
state). 



A multi-soliton is a solution of (NLS) built with solitons. More precisely, let A^ £ N \ {0,1}, 
uji,...,ujn > 0, 7i,...,7Ar e R, vi,...,vn e M'*, xi,...,xn G K'' and $i,...,$Ar £ H'^{R'') solutions 



of (|^) (with ijJq replaced by Wi, ..., wjy)- Set 



JV 



(3) Rj{t,x) -.^ R^-^u>j,fj,vj,xj{t,x), R{t,x) := } Rj(t,x) 



j=i 
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Due to the non-linearity, the function R is not a solution of (NLS) anymore. What we call multi 



soliton is a solution u of (NLS) defined on [To, +00) for some Tq £ R and such that 



lim \\u{t) -i?(i)||^i(Rd) = 0. 

i— f+00 ^ ^ 

In this paper, we are concerned with existence, non-uniqueness and instability of multi-solitons 
build on excited states, which we will refer to as excited multi-solitons. 

History and known results. Solitons and multi-solitons play a crucial role in understanding 
the dynamics of nonlinear dispersive evolution equations such as Korteweg-de Vries equations or 
nonlinear Schrodinger equations (see e.g. [p7| for a general overview) 



To fix ideas, consider the pure-power nonlinearity f{u) = juj^^^u. Equation (NLS) is L'^-critical 
(resp. subcritical, resp. supercritical) \i p — 1 + ^ (resp. p<l-f^,p>l-|-^). The soliton 
resolution conjecture states that, at least in the L^-subcritical case, a generic solution will eventually 
decompose into a sum of ground state solitons and a small radiative term, in some sense we will 
not try to make precise. However, this conjecture remains widely open, except when the equation 
is completely integrable (like the classical Korteweg-de Vries equation ut + Uxxx + uux = 0) and 
explicit solutions are known [g6| ^ . 

Nevertheless, multi-solitons based on ground states are supposed to be generic objects fo r larg e 
time; in contrast excited multi-solitons are believed to be singular objects of the flow of ( |NLS ). 



However, their existence shows that a global approach of the large time dynamics must take care 
of them. 

The first existence result of multi-solitons in a non-integrable setting was obtained by Merle 
[p2[ for multi-solitons composed of ground states or excited states for the L^-critical nonlinear 
Schrodinger equation. For multi-solitons composed only of ground states, the L^-subcritical case 



was treated by Martel and Merle [^ (see also Martel 29 for the generalized Korteweg-de Vries 
equation) and the L^-supercritical case by Cote, Martel and Merle [Q. No excited multi-solitons 
were ever constructed except in the L^-critical case and our result (Theorem ||) is the first in that 
direction: we construct excited multi-solitons based on excited states which move fast away from 
one another. 

Study of the dynamics around ground-states solitons and multi-solitons, in particular stability 
properties, has attracted a lot of attention since the beginning of the 80's (see e.g. [2], g, ^ ^, 
BM p% Eoi). The main result states that ground-states solitons are orbitally stable only in the 
Z?-subcritical case. 

So far, little is known about the stability of excited state solitons. All excited states are 
conjectured to be unstable, regardless of any assumption on the nonlinearity. For results on 
instability with a supercritical nonlinearity, see Grillakis [19| and Jones pJ in the case of real 
and radial excited states and Mizumachi for vortices [pSl |34[. Partial results in the L^-subcritical 
case are available in the works of Chang, Gustafson, Nakanishi and Tsai [^, Grillakis |g^ and 
Mizumachi [ pS] . 

Here we show that under a very natural assumption of instability of the linearized flow around 
one excited state, the excited multi-soliton is not unique, and unstable in a strong sense. 

Statement of the results. We make the following assumptions on the nonlinearity (recall that 

/(z)=5(|z|2)zforzGC). 

(Al) geC°([0,-Koo),K)nCH(0,+oo),R), g{0) = and lim^^o s.9'(s) ==0. 

(A2) There exist C > and 1 < p < 1 -I- ^ ii d ^ 3, 1 < p < +00 if d = 1,2 such that 

\s^g'{s^)\ s; CsP~^ fors^ 1. 
(A3) There exists Sq > such that F{sq) > ^. 

Remark 1. A typical example of a non-linearity satisfying (A1)-(A3) is given by the power type 
non-linearity f{z) = \z\p~^z with 1 < p < 1 -|- -j:^ if d ^ 3, 1 < p < -foo if d = 1, 2. 

Assumptions (Al)-(A3) guarantee that, except in dimension d ~ 1 where all bound states are 
ground states, there exist ground states and infinitely many excited states (see e.g. [g|, ^, |[ po|, |25|]). 
In particular, excited states can have arbitrarily large energy and L°°(R'')-norm. Note that every 
solution of (|l|) is exponentially decaying (see e.g. [Q). More precisely, for all $0 solution to (||) we 
have e^^l^l(|$o| + |V$o|) e L°°(R'*) for ah w < wq. 
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Assumptions (Al)-(A2) ensure well-posedness in H^{R'^) of ( NLS| ) , see e.g. (the equation is 



then iJ-'^-subcritical). In particular, for any uq G H^{M.'^) there exists a unique maximal solution 
u such that energy, mass and momentum are conserved. Recall that energy, mass and momentum 
are defined in the following way. 

1 r 

F{u)dx, 



E{u) 


— 2 ^"^ i2(R<i) ~ 


M{u) 


- IU/I|2 

— ll"llL2(Kd), 


P(u) 


= J?^^ / uS/udx 



Notice that (A3) makes the equation focusing. 

Our first result is the existence of multi-solitons composed of excited states as soon as the relative 
speeds Vj — Vk of the solitons are sufficiently large. 

Theorem 1. Assume (Al)-(A3). Let N eN\ {0,1}, and for j = l,...,N take uj.j > 0, jj G R, 
Vj € M"^, Xj € M'^ and $j e iJ^(IR'^) a solution of (|l|) (with luq replaced by loj). Set 

Rj{t,x) = i?$,,^„^^-,„,,,^-(t,x) := $,(x - Vjt - a;,)e^(^^ — 3l'^^l'*+-^*+^^). 
Let LOi, and u* be given by 

^* '■= ^^i^i'^J'-? "" l7---,Af}, t'* := q™^"'fl"j ~ Vkl'J^k = l,...,N,j ^ k} . 

Also introduce a :— sin I ^j,.^, ) (this constant appears naturally in Claim Ua). 
There exists wj :— Wjj($i, ..., ^n) > such that if v^, > a^^Vf^ then the following holds. 



There exist Tq E M. and a solution of (NLS) u E C([To, +oo), H^{M.'^)) such that for all t G [Tq, +cx)) 
we have 

N 1 

|ly(t)-^i?,(i)|Ui(K.)^e— *-'. 

We now turn to the non-uniqueness and instability of a multi-soliton. 

Assume that the flow around one of the Rj is linearly unstable, i.e. has an eigenvalue off the 
imaginary axis. As the Rj all play the same role, we can assume it is i?i. 

(A4) L = —iA + iuJi — idf{^i) has an eigenvalue A G C with p :— Me{X) > 0. 
This assumption is very natural if one expects Ri to be unstable. Actually, (A4) holds for any real 
radial bound state in the L^-supercritical case (see [|9J). For excited states, (A4) is believed to 
hold for a wide class of non-linearities. 



Under assumption (A4), we are able to construct a one parameter family of solutions to (NLS) 
that converge to the soliton Ri as time goes to infinity, as described in the following Theorem. 

Theorem 2. Take ioi > 0, 71 G M, ui G R'^, Xi G R'^ and $1 G H^{R'^) a solution of (0) (with ujq 
replaced by loi). Set 

i?i(i,x)=i?*„^„^„,„„,,(i,x):=$i(a;-«ii-a;i)e^(^-^-^l''il'*+'^i*+^i). 

Assume g is C°° and (Al)-(A4) are satisfied. 

There exists a function Y{t) .such that \\Y{t)\\ffi/-^d\ ^ Ce^f* and e''*||y(i)||^i is non-zero and 
periodic (here p is given by (A4-) and Y(t) is actually a solution to the linearized flow around Ri, 



see (|26|), (27) j. For all a G M, there exist Tq G M large enough, a solution Ua to ( NLS ) defined on 
[To, +00), and a constant C > such that 

yt ^ To, \\uait) - i?i(t) - ay(t)||^i(„d) < Ce-^P'. 

In particular. Theorem implies that the soliton Ri is orbitally unstable, as precised in the 
following corollary. 

Corollary 2. Under the hypotheses of Theorem 0, i?i is orbitally unstable in the following sense 



Let CT > 0. There exist e > 0, (r„) C (0,+oo), (uo,„) C i/^(R'*) and solutions (u„) of ( [NLS| ) 
defined on [0,T„] with u„(0) = Uo,n such that 



lim ||uo,„ - i?i(0)||/f,.(Rd) = and inf ||u„(r„) - e* $i(- - y)||L2/Rd) > e /or a^/ n G N. 
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From Theorem we infer the existence of a one parameter family of muhi-sohtons. As a 
corohary, we obtain non-uniqueness and instabihty for high relative speeds multi-solitons. 

Theorem 3. Let N eN\ {0, 1}, and for j = 1, ..., N take lUj > 0, 7^ G R, Vj G R'^, Xj G R'^ and 
$j e H^{R'^) a solution of ffl) (with ujq replaced by Wj). Set 

Let Ui := g min {\vj — Vk\;j,k — 1, ..., N,j 7^ k} . Assume g is C°° and (Al)-(A4) are satisfied. 

There exists Vi^ := wt]($i, . . . , ^n) > such that if u* > Vi^ then the following holds. 
There exists a function Y(t) such that \\Y{t)\\jjit]^d\ ^ Ce"''* and e''*||l"(i)||/fi is non-zero and 
periodic (here p is given by (A4-) and Y{t) is actually a solution to the linearized flow around Ri, 



see (Pq), (|27|)j. For all a €R, there exist Tq (zR large enough, a solution Ua to (NLS) defined on 
[To, +00), and a constant C > such that 

N 

yt^To, \\ua{t) -Y^Rjit) - aY{t)\\H^M-^) ^ Ce-^"- 

Remark 3. Notice that, in Theorem [J, if for a, 6 G M we have a ^b, then Ua^ u^,. Indeed, for t 
large enough we have 

\\ua{t) - WfeWllHMR^) ^ |a - b\\\YmHHR^) ~ 2Ce~^P'. 
Since e''*||y(t)||//i is non-zero and periodic, this implies that Ua ^ Ub if a, ^ b. 

Corollary 4. Under the hypotheses of Theorem [J, the following instability property holds. Let 
a ^ 0, there exists e > 0, such that for all n Cz M\ {0} and for allT (^R the following holds. There 



exists In, Jn G M, r ^ /„ < Jn and a solution Wn G '^([/n, Jn],II^{R'^)) to (NLS) such that 

N 

lim ||u;„(/„)-i?(/„)||H.(Rd) =0, and inf ||w;„( J„)- V $j(a;-yj)e'(3«^-^+'»j)||^2.jj^. ^ e. 

] = 1,...,N 

Remark 5. The fact that instability holds backward in time (i.e. with Jn < In) is an easy 
consequence of Theorem H. Hence the difficulty in Corollary is to prove instability forward in 
time. 

Remark 6. The classification of multi-solitons is now complete for the generalized Korteweg-de 



Vries equations (see |12, 29 and the references therein). In particular, uniqueness holds in the 
subcritical and critical cases, whereas in the supercritical case the set of multi-solitons consists in a 
N -parameters family. To the authors knowledge, no uniqueness nor classification result is available 
yet for multi-solitons of nonlinear Schrodinger equations. 

Scheme of proofs and comments. Our strategy for the proof of the existence result (Theorem n|) 
is inspired from the works O, ^ |32| : we take a sequence of time T„ — >■ +00 and a set of final 



data Un{Tn) = R{Tn). Our goal is to prove that the solutions Un to (NLS) backwards in time 
(which approximate a multi-soliton) exist up to some time Tq independent of n, and enjoy uniform 
H^{R'^) decay estimates on [To,T„]. A compactness argument then shows that (u„) converges to 
a multi-soliton solution to (NLS|) defined on [Tq, -l-oo). 



As in p3l pOl, the uniform backward iJ^(R^)-estimates rely on slow variation of localized 
conservation laws as well as coercivity of the Hessian of the action around each component of 
the multi-soliton. However, this Hessian has negative "bad directions" on which it is not coercive. 
When dealing with ground states, these were ruled out either by modulation and conservation 
of the mass (as in |Q) or with the help of explicit knowledge of eigenfunctions of the operator 



corresponding to the linearization of ( NLS ) around a soliton (as in ||13|). In both cases, this could 
be done only because of the knowledge of precise spectral properties for ground states; this does 
no longer hold when dealing with the more general case of excited states. 

Our remark is that the Hessian fails to be iJ-'^(K.'^)-coercive only up to a L^(M'^)-scalar product 
with the bad directions. Hence the first step in our analysis is to find uniform L^(IR'')-backward 
estimates without the help of the Hessian. This rules out the "bad directions" and we can now take 
advantage of the coercivity of the Hessian to obtain the _ff^(R'')-estimates. The main drawback of 
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our approach is that the bootstrap of the L^(M'^)-estimates requhes that the soUton components 
are well-separated. Thus we have to work with high-speed solitons. 

To obtain the one parameter family of Theorem 0, we rely on a fixed point argument for smooth 
functions exponentially convergent (in time). This is possible because we now assume smoothness 
on the non-linearity. The main difficulty is to construct a very good approximate solution to 
the multi-soliton. Actually we build such a profile at arbitrary exponential order. This method is 
inspired by ||l^, ll^, EtI Esj in the case of a single ground state, for the nonlinear wave or Schrodinger 
equations. It was also recently developed by Combet ||ll|, ^ for multi-solitons in the context of 
the L^-supercritical generalized Korteweg-de Vries equation. 

However, an important difference in our case is that we consider excited states, and the linearized 
flow around them is much less understood than that around a ground state soliton. For example, 
to our knowledge, the exponential decay of eigenfunctions was not known in general (see ||23| for a 
partial result). We prove it inQ see Proposition^ Also, the unstable eigenvalue has no reason to 
be real, and this will make the construction of the profile much more intricate than in the ground 
state soliton case. This is the purpose of Proposition |2[ Once the approximation profile is derived, 
the proofs of Theorem and follow from a fixed point argument around the profile. 

The paper is organized as follows. In Section 2 we prove Theorem |l|. Section 3 is devoted to 
the proofs of Theorems and p. In H we prove the exponential decay of eigenfunctions for matrix 
Schrodinger operators and in p we prove Corollaries and H. 

Acknowledgement. The authors are grateful to the unknown referee for valuable comments and 
suggestions. 

2. Existence 

In this section, we assume (Al)-(A3) and suppose we are given N E N\{0, 1}, and for j — I, ..., N, 
ojj > 0, 7j G M, Vj e K'', Xj £ W^ and $j e H^{R'^) a solution of (Q) (with cjq replaced by ujj). 
Recall that 



i?^-(t,a;) = $j-(x-Wji-Xj-)e'(5"^-^-|l''.l'*+'^.t+7.)^ 
-mm{wj,j = 1,...,7V}, «* = g 



^* = o^'^i^i'-? "= l,---,^}, «* = nmin{|uj -~Vk\;j,k = l,...,N,j ^ k} 



and a := sm y j^^^^l^ 

2.1. Approximate solutions and convergence tovirard a multi-soliton. L et (T„ )n:gi C M 
be an increasing sequence of time such that T!„ — > -|-oo and (u„) be solutions to ( NLS| ) such that 
Un{Tn) = R{Tn). Wc Call u„ an approximate multi-soliton. 
The proof of Theorem |l| relies on the following proposition. 

Proposition 7 (Uniform estimates). There exists wj := wj($i, ..., <I>jv) > such that ifv^, > a~^Vf^ 
then the following holds. There exist uq € N, Tq > such that for all n ^ uq every approximate 
multi-soliton u„ is defined on [To,T„] and for all t e [To,r„] we have 

(4) ||zi„(t)-i?(t)||Hi(R.)^e-"-*'^**. 

In this section, assuming Proposition R, we prove Theorem 1 by establishing the convergence of 
the approximate multi-solitons m„ to a multi-soliton u existing on [Tg, -f-oo). Our proof follows the 
same line as in [^ pO[ . 



From now on and in the rest of section 2^ we assume that w* > a ^wj, where wj is given by 
Proposition 0. 



Since the approximate multi-solitons Un are constructed by solving ( NLS ) backward in time, to 
prove Theorem |l| we first need to make sure that the initial data Un{To) converge to some initial 
datum uq. 

Lemma 8. There exists uq G H^lM."^) such that, possibly for a subsequence only, u„(To) — >■ uq 
strongly in iJ'*(M'^) as n — > +oo for any s G [0, 1). 

Lemma H is a consequence of the following claim. 
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Claim 9 (L^(M'*)-compactncss). Take S > 0. There exists rs > such that for all n large enough 
we have 

(5) f \uniTo)\^dxi^S. 

J\x\>rs 

Proof. Let n be large enough so that the conclusions of Proposition ^ hold. Let Ts be such that 
^-au^v.Ts ^ /| rpj^gj^^ ^y Proposition 0, we have 

(6) \\uniTs) - R{Ts)\\h^r'') =^ Y 4' 
Let ps be such that 

(7) / \R{Ts)\'dx < ^. 



\x\>ps 



From (||)-(0) we infer 

(8) 



/ \Ur,{Ts)\^dx < ^. 



We define a C^ cut-off function r : M ^ R such that t{s) = if s ^ 0, t{s) = 1 if s ^ 1, t{s) e [0, 1] 
and |t'(s)| ^ 2 if s G [0, 1]. Let ks to be determined later and consider 

T(t) := / |«„(t)|V ("M^^'j dx. 

To obtain (|h from M) we need to establish a link between T(To) and T{Ts). Differentiating in 
time, we obtain after simple calculations (see e.g. |^, Claim 2]) 

T (t) — Jl^ / UnVUn ■ T-rT d.x . 

US jRd \x\ \ Kg J 

Since |lMn(^)|lHi(R<i) is bounded independently of n and i, there exists 

Co := sup sup ||w„(i)||^i(jjd) > 

TiGNte[To,T„] 



such that 



|T'(t)| «; ^. 

Kg 



Choose Kg such that ^§^Ts < |. Then, by integrating between Tq and T^ we obtain 

(9) T(To)-T(r,)sC^. 

From (ph we infer that 



T(T5) = / \uniTs)\'T L^^ll \dx^ \uniTs)\'dx < - 

JWL'^ \ ^^S / J\x\>ps ^ 

Combining with (M) we obtain 

T(To) ^ 5. 
Now set rs :— Kg + ps- Then from the definition of r it is easy to see that 

\ujTo)\^dx^TiT„)^S, 

\x\>rs 

which proves the claim. D 

Proof of Lemma ^. Since u„(ro) is bounded in H^(R'^), there exists uq G H^(R'^) such that up 
to a subsequence u„(Tb) -^ uo weakly in H^{M.'^'). Hence, u„(To) — i' wo strongly in Lf^^{M.'^) and 
actually strongly in L^(IR'*) by Claim 0. By interpolation wc get the desired conclusion. D 
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Proof of Theorem \^. Let uq be give n by Lemma | and let u S C{[To,T*), H'^{R'')) be the 
corresponding maximal solution of ( NLS| ). By (Al)-(A2), there exists < ct < 1 such that 



1< P < 1 + ra and 

|/(zi) - f{z2)\ ^ C{1 + |zi|P-i + \z2r^)\zi - Z2\ for all zi,Z2 € C. 



This implies that the Cauchy problem for (NLS) is well-posed in H"{W^) (see R S])- Combined with 
Lemma this implies that u„(i) — >■ u{t) strongly in H'^{W^) for any t G \Tq^T*). By boundedness of 
u„(i) iniJi(E'^), we also have M„(i) -^ u{t) weakly 'yd. H^ {W^) for anyi G [T^.T*). By Proposition 0, 
for any t £ \Tq^T*) we have 



(10) \\u{t) - i?(t)||Hi(Rd) ^ lini|nf ||u„(i) - i?(i)||^i(Rd) ^ e 



-acj^ -u^i 



In particular, since R{t) is bounded in H^{W'-) there exists C > such that for any t G [Tq, T""*) we 
have 

(11) \\u{t)\\m(K^) ^ e-""?-* + ||i?(i)||Hi(M^) < C. 

Recall that, by the blow up alternative (see e.g. 0), either T* = +oo or T* < +oo and 
limi_j.T*||w||//i(Rd) = +00. Therefore ( pT| ) implies that T* = +cx). From ( pl]| ) we infer that for 
all t E [Tq, +oo) wc have 

h(t)-i?(0||ffi(R.)^e-"-*''**. 
This concludes the proof. D 

2.2. Uniform backward estimates. This section is devoted to the proof of Proposition]^. This 
proof relies on a bootstrap argument. Indeed, from the definition of the final datum w„(T'„) and 
continuity of u„ in time, it follows that (||) holds on an interval [i^,T„] for f^ close enough to r„. 
Then the following Proposition nO shows that we can actually improve to a better estimate, hence 
leaving enough room to extend the interval on which the original estimate holds. 

Proposition 10. There exists wj := Wjj($i, ..., ^n) > such that if Vi, > a^^wj then the following 
holds. There exist ng G N, Tq > such that for all n ^ rip every approximate multi-soliton u„ is 
defined on [To,T„]. Let t^ G [To,T„] and n ^ uq. If for all t G [t^,Tn] we have 

(12) ||u„(i)-i?(i)||^i(R.)<e-""*'^*' 
then for all i G [t^ , Tn] we have 

(13) h„(i)-i?(i)||^i(R.)S^ie-"'^*''**. 

Before proving Proposition |lO|, we indicate precisely how it is used to obtain Proposition [^. 

Proof of Proposition |^. Let Tq, uq and wj be given by Proposition ^ assume v^ > a^^v^, and let 
n ^ no. Since u„(T„) = R[Tn) and u„ is continuous in H^{W^), for t close enough to T„ we have 

(14) |li.„(t)-i?(t)||Hi(R.)S=:e-""*'"**. 
Let t^ be the minimal time such that dlj) holds: 

<^ := min{r G [ro,T„]; (H) holds for aU t G [T,Tr,]}. 
We prove by contradiction that t^ = Tq. Indeed, assume that t^ > Tq. Then 

h„(it)-i?(ft)||^,^^,^^e-""*''**' 
and by Proposition |lO| wc can improve this estimate in 

ii«„(it)_i?(it)||^,^^,^^i -""*^**^ 



Hence, by continuity of u„(t) in H^{W^), there exists Tq ^ t^ < t^ such that (14) holds for all 



i G [i"'',i^]. This contradicts the minimality of t^ and finishes the proof. D 
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The proof of Proposition M is done in two steps. First, assuming (Q2|) we prove that we can 
control the L^(K'^)-norm of (m„ — R). To obtain the full control on the i?^(R'^)-norni of (un — R) 
as in (113) we use the linearization of an action-like functional. This linearization is coercive (i.e. 
controls the H^ (R'')-norm) up to a finite number of non-positive directions that can all be controlled 
due to the L^(]R'')-estimate. 

Let To > large enough and fix n € N such that T„ > Tq. For notational convenience, the 
dependency on n is understood for u and we drop the subscript n. Set v := u — R. Let t^ G [Tq, T„] 
and assume that for all t G [t^ ,Tn] we have 

||i;(i)||ffi(M.)^e-"-*^**. 
Step 1: L'^{W^)-control. 

Lemma 11. For all K > and ?7i e N \ {0} there exists uj = v^(K,m,^i, ...,^n) > such that 
if Vi, > a^^v^ then for all t £ [t^jTn] we have 

Notice that the reason why we introduce such K and m will appear later in the proof. 

Proof. First note that by identifying C to R^ and viewing / : R^ — ;> R^ we can consider 

df{z).w = g{\z\^)w + 2^4ziu)g'{\z\^)z. 

The function v satisfies 

ivt + Cv + M{v) = 0, 
where 

Cv := Av + df{R).v 
and the remaining nonlinear term Af{v) verifies 

I (i7V(w),w)i2(Rd) I s; v{\\v\\Hiim^))\\v\\H^jidy 

where 77 is a decreasing function satisfying 77(5) — ^ when s — > 0. Take any t G [v ,Tn]. We have 
1 d „ „, 



We have 



2Jt'"'^'( 



{iCv,v)^2 



{vuv)^2(^a) = (i£u,u)^2(Rd) -f (iA/'(t;),w)i2(Rd) 



i{l^v + df{R).v))vdx, 
i {Av + gi\R\^)v + 2g'{\R\^)!^e{Rv)R) vdx, 
S^. I ii-\Vv\^ + g{\R\^)\v\^ + 2g'{\R\^)^4Rv)Rv)dx, 
2g' {\R\^)^e{Rv)j^^{Rv)dx. 



Therefore, 

|(z/:«,7;)^2(K.)K / 2\g'{\R\')\\R\'\v\^dx. 

where this last constant C'c depends only on g and ||i?||^oo(jjd-). By the bootstrap assumption on 
V, this implies 

In addition, it is easy to see that 



(i7V(u),u)^2(Rd) I < viMm{Rd))\\v\\H^Tg.d^ < r]{- 



e 
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In short, if Tq is large enough so that 77(6"""* "**) ^ ^, we have obtained that 

d_ 
di' 

Therefore, by integration between t and T„ we get 



nlkllL^fRd) 



^2Cce 



— 2auj^ v^t 



(15) Mml.^^.) - ||«(T„)|li.(K.) ^ -^(e-^-?"** _ e-2a^?..T„ 

Now, we take Vf^ such that 






If Wi > a -^wj and since f (T„) = we get from ( fi5| ) that 

1 i 

which is the desired conclusion. D 

Step 2: H'^iW^)- control. The idea of the second step of the proof of Proposition |l^ is reminiscent 
of the technique used to prove stability for a single soliton in the subcritical case (see e.g. 
pi| , p4 p7| , p9| , Eoi). Indeed, it is well-known that the linearization of the action functional Sq 
(see the definition of Sq p. 0), whose critical points are the solutions of dl^), is coercive on a 
subspace of H^{W^) of finite codimension in i^(R'^). At large time, the components of the multi- 
soliton are well-separated and thus it is possible to localize the analysis around each soliton to 
gain an if ^(R'^)-local control, up to a space of finite dimension in L'^iM.'^). But due to Lemma O 
we are able to control the remaining L^(R'^)-directions, hence to close the proof. The idea of 
looking at localized versions of the invariants of (NLa) was introduced in [pa and later developed 



in [tl3l E3, M, pll]. We shall therefore be sketchy in the proofs, highlighting only the main differences 
with the previous works. 

We start with the case of a single soliton. 

Lemma 12 (Coercivity for a soliton). Let wq > 0, 70 G R, a^Oi vq G R'' and a solution $0 G H^{W'') 
of (I). Then there exist Kq = iCo(*o) > 0, i^q E N \ {0} and X^, ...,Xo° £ L'^{R'^) such that for 
k — 1, ..., z^o we have \\XQ\\i^2md\ = 1 and for any w G H^{M.'^) we have 

Ikllffi(R^) ^ KoHoit,w) + Koy2{w,X^{t))l forallteR, 






wher 



Xl;{t) ■-e''^ho-^-i\vo\'t+u:ot+^o)xk(^^ _ ^^^ _ 3.^)^ 



Hoit,w) :=||Vw||^2(Rd) + [ujo + ^— ) ||w||;^2(Rd, - vo ■ J^^^ / wVwdx 

2\L,,|2 I o„'/lD |2\,52 /D „T-.^2^ 



{g{\Ro\^)\w\^ + 2g\\Ro\^)^4Rowy) dx, 

and Ro{t,x) is the soliton given by (0). 

Lemma |l^ follows from standard arguments. We included a proof in y for the reader's 
convenience. 

We introduce now the localization procedure around each component of the multi-soliton. 
We begin by the selection of a particular direction of propagation. 

Claim 13. Let < a < sin I ^ \\r('i\ ) ■ ^^s*^ there exists an orthonormal basis (ei, ..., Cd) of 
R"* such that for all j,k = 1, ..., N , we have 

\{vj -Vk,ei)Rd\ ^ a\vj - Vk\. 
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Proof. For j ^ k, set Vjk ■— \y'._y''\ ■ The claim will be proved if we show that the measure of the 

set 



A:= U {weS''-\\{v,k,wM^a} 



j,fc=l,.-.,iV 

is smaller than the measure of the surface of the unit sphere S'^~ ^ . 

Take j,k = 1,...,N; j 7^ k. Without loss of generality, assume that Vjk = (1,0,..., 0). Take 
w G S'^^^ and let {9i, ...,9d-i) be the spherical coordinates of w. Then we have 

{vjk,w)s,d = cos 6*1. 

Therefore, after easy calculations we get 



^({w e S , \{vjk,w)jgd\ ^ a}) < 2arcsin(a) 



d-l 

n 2 



d-l 

^ and Zd-i; is the area 
r(V-) 

subadditivity of the measure this leads to 



where fi is the Lebesgues measure on S and J',/-!, is the area of the {d — 2)-unit sphere. By 



d-l 

TT 



Now, remember that 



This implies 



/.(A) s; N{N - 1) arcsin(a)— ^— . 



0<a<sin' ^ ^ ^ ' 



iV(A^-i)r(f) 



d-l 
TT 2 



/i(A) ^ 7V(iV - 1) arcsin(«)-^^ < Wd\ = ^<^ )■ 
Therefore /x(§''^^ \ A) > and we can pick up ei e S"^^^ such that for all j, k — 1, ..., N, we have 

Completing ei into an orthonormal basis (ei, ..., e^) of M.'^ finishes the proof. D 



By invariance of (NLS|) with respect to orthonormal transformations we can assume without 



loss of generality that the basis (ei, ..., e^) is the canonical basis of M.'^. Up to a changes of indices, 
we can also assume that vl < ... < v^^ where the exponent 1 in v^ denote the first coordinate of 
Vj ^{v],...,v'j). 

Let V) : K ^ M be a C°° cut-off function such that tp{s) = for s < -1, V(s) e [0, 1] if s e [-1, 1] 
and ip{s) — 1 for s > 1. We define 

m, := -{v]_^ + v]) for j = 2, ..., N, 

ipiit,x) := 1, Vj(^,a;) := V(^(a;^ ^ mjt)) for j = 2,...,N. 

Then wc can define 

We introduce localized versions of the energy, charge and momentum. For j = 1, ..., TV we define 

Ejit,w):^- \Vwf<j)jdx- F{w)(t)jdx, 



Mj{t,w) -.^ / \w\4)jdx, Pj{t,w) -.^ .y^-ri. I {Vw)w<j)jdx. 

JwLd jRd 

We denote by Sj a localized action defined for w G H^(R'^) by 

Sj{t,w) := E,{t,w) + i (^LJ, + ^^ M,it,w) - \, ■ P,{t,w) 
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and by Hj a localized linearized defined for w E H^{M.'^) by 

Hj{t,w):^ / |Vwp0jda;- / {g{\Rj\'^)\w\^ + 2g' {\Rj\'^)^^{Rjwf) (fjjdx 



OJj H ^ — 1 / \w\'^(j)jdx ~ Vj ■ ,_^>rt \ w\/w(j)jdx. 

We define an action-like functional for multi-solitons 

N 

S{t,w) ■.= ^Sj{t,w) 
i=i 
and a corresponding linearized 

N 

n{t,w):=Y,Hj{t,w). 

We have the following coercivity property on H. 

Lemma 14 (Coercivity for the multi-soliton). There exists K — K{^i, ...,^n) > such that for 
all t large enough and for all w £ H^{M.'^') we have 

N t^i 

J=l 1=1 



li 



where {vj), {Xy) are given for each Rj by Lemma 

Proof. It is a consequence of Lemma E2| (see [pll Lemma 4.1]). D 

Lemma 15. The following equality holds 

5,(t,u(i)) = 5,(i,i?,) + iJ,(i,T;) + 0(e-3"-?"**) + o(||z;||^,(„.)). 
The proof relies on the following claim. 
Claim 16. For all x E M'' and j, k — 1, ..., N the following inequalities holds. 

{\Rk{t,x)\ + \VRk{t,x)\)cj,,it,x) ^ Ce-2"-* ''**e-^l^-'"=*-^^l for j ^ k, 

i\R,it,x)\ + \VR,it,xMl ^ q^,it,x)) ^Ce-^^-i^^^'e-^l"-"^*-"^!. 

Proof. The claim follows immediately from the support properties of 0^ , the definitions of io^, and 
v^, and exponential decay of $ j . D 

Proof of Lemma \1^. The proof is done by writing u{t) = R{t)+v{t) and expanding in the definition 
of Sj. We start with the terms of order in v. By Claim [iq we have 

(16) Sj{t,R) ^ Sj{t,Rj)+0{e~^°"^* ""*'). 

We now look at the terms of order 1 in w. Still by Claim PlgL taking in addition into account that 
1 ' ' — ' 

||w||^i™d) = 0(e~""* "**) and remembering the equation solved by Rj (see (43)) we obtain, 



1-^11 y- \ /—\ / — "irv/.i -^ 1 



(17) (5j(t,i?),«) = (5j(t,i?,),«)+0(e-3"-*"**)=0(e-3"-* ''**), 

(18) (5;(i,i?)t;,t;) ^iJ,(i,t;) + 0(e-3"-*"-')+o(||z;||?,i(„.)). 
Gathering (|l^)-(|l^) we obtain the following expansion 

S,{t,uit)) = 5,(i,i?,) + iJ,(i,t;) +0(e-3""?"*') + o(||t;||^,(j,.)), 
which concludes the proof. D 

We can now write a Taylor-like expansion for S. 
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Lemma 17. We have 

S{t,u)~S{t,R)=n{t,v)+o{\\v\\l^^^,^) + 0{e-^"'^*''*') 
Proof. In view of Lemma n^ all we need to prove is 

N 1 

5(t,i?)=^^,(t,i?,)+0(e-3"-? ''**), 
which follows immediately from Claim hq. 



n 



Lemma 18. The following estimate holds. 

dS{t,u{t)) 



dt 



^ ~2au:'^Vtt 



^^e 



Proof. We remark that 

S{t,w)=E{w) + J2(l 



N 



J = l 



^j+^-^)M,{t,w)-^v,-P,{t,w) 



Since the energy E is conserved by the flow of ( NLS| ), to estimate the variations of S{t,u{t)) we 
only have to study the variations of the localized masses Mj(t,u{t)) and momentums Pj{t,u(t)). 
Take any j — 2, ..., N. We have 



1 

^7t 

Define Ij :— [rrijt — ^/t, lUjt + ^/t] x 

d 



2X^ + mjt\ ,, , 1 



At 



ip'{—^{x^ — mjt))dx. 
\t 



Q^ I ,u(i)|Vj(i,a:)(ix 



■^ ^ . From (ny) and the support properties of ip we obtain 
C 



€ -^ / \Vu\^ + \u\^dx. 
yt Ji- 



Similarly, for the first component of Pj we have 

I d [ _ 

(20) o^ / udiuipjdx = 

I at j^d 



Vt L ('^''''' " ^C"'')''''' + ^("'^ ~ "^1"^-^) ^'^^y - "^^"^)) 

-^^\u\^r'{^y~m,t)dz 
Combining (EOh with the support properties of -0 and (A1)-(A2) we obtain 



— / udiwpjdx 



s;-^ I / \Vu\^ + \u\^dx 



\\7u\^ + \u\^dx 



Similar arguments lead for fc ^ 2 to 

d 



dt 



udkuipjdx 



< 



C 



ViJi, 



iVwP + luPdx. 



Now, we remark that 



/ {\Vu\^ + \u\^)dx^ f |Vi?|2 + |i?|2da;+||u-i?||^i(K<*) 
Recall that by hypothesis we have 

In addition, the decay properties of each ^k and the definition of Ij imply 

1 



/ (|Vu|2 + \uf) dx ^ Ce-2""* 
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Consequently, 



- / \uit)\'^,{t,x)dx 



dt ' ^ 



^ —pe * ■ 



Note that the previous inequahty is trivial for i ~ \ since -ipi ~ 1 and the mass and momentum 
are conserved. Plugging the previous into the expressions of Mj and Pj gives 



iM,{t,u)+P,it,u)) 



C 



-2qlj/ v^t 



and the desired conclusion readily follows. 



D 



Proof of Proposition u3. Let K = K{^i^ ..., ^n) and m := X]7=i ^j t)e given by Lemma nj. Since 
||XJ^(t)||^2(Rd) = 1 for any t,j,k, by Lemma |ll|, there exists V(^ = Vf^{^i, ...,^n) such that if 
w* > a^^fjj we have for j — 1, ..., N, k — 1, ..., Vj that 



(21) 



(i;(i),4W)^.(K.)^lkWlli.(M^) 



€ 



1 



2mK 



Using Lemma n^ we obtain 



(22) 



S{t,u{t))-S{Tn,u{T„))^ 



dS{s,u{s)) 



ds 



-2aij^ v^t 



c 



ds < ^e-^""* ''**. 



Vi 



Note that since u„(T„) — i?(r„) we have 

(23) 5(T„,7/(T„))-5(T„,i?(r„)) = 

By Lemma0,(|2|)-(|2|) imply 



(24) 



nt,v) ^ — ^1 — +o(ii«ii^i(R.)). 



Combining (pi|)-(p4|) and Lemma |lj we get 

and we easily obtain the desired conclusion if To is chosen large enough. 



n 



3. Non-uniqueness and instability 

In this section, we assume g e C°° and (Al)-(A4) are satisfied. We take N E N WO, 1}, and 
for j = l,...,iV, ujj > 0, 7j e M, Vj e M'', x^ e 
replaced by ujj). Recall that 



and $j e i/i(M'') a solution of (m) (with wq 



w* = 7:min{wj,j = 1,...,A^}, w^ = -min{|t;j -'yfel;^,^ = l,...,N,j ^ fc} 



3.1. Construction of approximation profiles. Since ( |NLS ) is Galilean invariant, we can 
assume without loss of generality that vi = 0, 71 — 0, xi — 0. For notational brevity we drop 
in this subsection the subscript 1 indicating that we work we the first excited state. Hence we will 
write (in this subsection only) Ri{t,x) = R{t,x), $1 — <&, etc. 

Note first df{z).w = g{\z\'^)w + 23ie{zw)g' {\z\'^)z is not C-linear. This is why we shall identify C 

with R'^ and use the notation a + ib = I , ] (a, & e M), so as to consider operators with real entries. 

Given a vector w G C^, we denote w+ and v^ its components (so that if v represents a complex 
number, 1;+ is the real part and v~ the imaginary part). To avoid confusion, we will denote with 
an index whether we consider the operator with C, M^, or C^-valued functions. 
Thus, as we consider 



^cv = -iAv-idf{R).v, 



-iAv + iujv — idf{^).v, 
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and the non-linear operators 

^c{v) = if{R + v)- if{R) - idfiR).v, 

^c{v) = e-'"^'yy{e"^'v) = i/($ + v) - i/($) - idf{^).v, 

then for instance 

^""'[v- -[-A + LO-I+ J [v- 



with $+ and $ the real and imaginary parts of $ and 

J = 2$+$-5'(i$n, /± = g(i$|2) + 2$± V(i*n- 

Now Lk2 is as an (unbounded) R-linear operator on iJ^(R'',R^) -^ L^(M'',M^). So as to have some 
eigenfunctions, we can complexify, and we are interested in Lc2 : H'^{M.'^, C^) -^ L'^{M.'^, C^), which 
is a C-linear operator with real entries. 

Let a > be the decay rate given by Proposition ^ for eigenfunctions of L with eigenvalue A 
(see (A4)). Possibly taking a smaller value of a, we can assume a £ (0, a/cj)- For K = M, M^, C or 
C^, denote 

(25) ^(K) = {v e H°°{R'^,K)\ e"l^l|i:i"w| G L°°(R'') for any muhi-index a}. 

We have gathered in the following proposition some properties of Lc2 that shall be needed for our 
analysis. 

Proposition 19 (Properties of Lc^)- 

(i) The eigenvalue X = p + i9 £ C of _Lc2 can be chosen with maximal real part. We denote 

Z{x) = „_\ i J G iJ^(R'^,C^) an associated eigenf unction. 

(ii) $ G ^(R2) and Z G ^^(C^). 

(iii) Let ^ ^ Sp(Lr2), and A G J^(C^). Then there exists a solution X G J^(C^) io 
(L — pI)X = A, and (L — pl)^^ is a continuous operator on Jif{C'^). 

Exponential decay of eigenvalues of L is a fact of independent interest. Hence we have stated the 
result under general assumptions in the |A| (see Proposition g5|). Notice that we treat all possible 
eigenvalues (in particular without assuming | jf^^Aj < cj, as it is the case for example in p3[). 

Proof, (i) It is well known that the spectrum of Lc2 is composed of essential spectrum on 
{*2/i y G K, \y\ ^ w} and eigenvalues symmetric with respect to the real and imaginary axes (see 
e.g. 1^, ^). The set of eigenvalues with positive real part is non-empty due to (A4). As L£2 is a 

compact perturbation of . n ] there exists an eigenvalue A with maximal real part. 

y — A -^ LJ \) J 

(ii) Exponential decay of $, V<i> is a well-known fact (see e.g. Q]). Then using the equation 
satisfied by $, one deduces that $ G Jf (R^). The decay and regularity of the eigenfunction Z rely 
essentially on the decay and regularity of $. Therefore, we leave the proof to ^, Proposition ^ 
and Proposition ^. 

(iii) Regularity of X follows from a simple bootstrap argument. For the exponential decay, we 
use the properties of fundamental solutions of Helmoltz equations (see Proposition ^0|). D 

To conclude with the notations, we define the decay class 0(x(i)), which we will use for functions 
decaying exponentially in time. 

Definition 20. Let S, G C°°(R+, i/°°(R'')) and x : R+ -^ (0, -t-oo). Then we denote 

£,{t) = 0{x{t)) as t -^ +0O, 
if, for all s ^ 0, there exists C(s) > such that 

Vf^O, \\m\\H^iR-)^C{s)x{t). 

Let Yi := ^.{Z) = (J[f-!) and Y2 :- J..{Z) = (J"K-}) ■ Then ¥^,¥2 G Jf{R^), and 

Lr2¥, = p¥,-9¥2, 
^2 = e¥i+p¥2. 
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Denote 

(26) Y{t) ^ e^P\coa{et)Yi + sm{et)Y2). 
Lemma 21. The function Y verifies for all t CzM. the following equation. 

(27) dtY + L^2Y = 0. 
Proof. Indeed, we compute 

dt{e-P\cos{9t)Yi+sm{9t)Y2)), 

= e-P^ {{~pcos{et) -9sm{9t))Yi + {- p sm{9t) + 9 cos{et))Y2) , 
LK2{e''"{cos{9t)Yi + sm{9t)Y2)), 

= e~''*(cos(6lt)iYi + sm{9t)LY2)), 

= e-P\cos{9t){pYi - 9Y2) + sin{9t){9Yi + pFa), 

= e-P\{pcos{9t)+9s\n{9t))Yi + {ps\n{9t) - 9 cos{9t))Y2 . 

So that {dt + LK2){Y{t)) =0. D 

Proposition 22. Let Nq eN andae M. Then there exists a profile W^" € C°°{[0, +00), Jf (M^))^ 
such that as t ^ +00, 

dtW^° + Lm2W^° = ^R2(W^^°) + 0(e-''(^°+i)*), 

and W^°{t) = ay(i) + 0(e-2p*). 

Remark 23. Notice that W^° {t, x) is a real valued vector. If we go back and consider W'^" as a 
function taking values in C, we then have, by definition of ^ , with U^°{t) = R{t) + e'"*VF^"(i), 

idtU^° + AU'^" + f{U^") = ©(e-P^^o+i)*). 

For the proof of Proposition E2l we write W for W^° (for simphcity in notation) and we look 
for W in the fohowing form 

No / k 

(28) W{t, x)^Yl e"''^* Y. ^J-^(^) cos(j0i) + B,^k{x) sm{j9t) 

fc=l \j=0 

where Aj k — ( J-' I and Bj k = ( j^-' I are some functions of .J^i^) to be determined. 

We start by the expansion of ^{W). 
Claim 24. We have 

No K 

J^^2{W) = J2 e"""* E {^jA^) cos(j0i) + B,, „(a;) sin(j0i)j + ©(e-'^^+i)"*) 

where Aj^^,, _Bj ^ € j?^(M^) depend on Ai^n and Bin only for I ^ n ^ k — 1. 

Proof. First we use a Taylor expansion. Due to smoothness of / and $ E ^(M^), and as ^r2 is 
at least quadratic in v, there exists a polynomial P/Vg G ^(R^)[X, F] with coefficients in Jf{M.'^), 
and valuation at least 2, such that : 

where Pj_„, Qj^„i e ^(R). 

Consider now the term W^W!"~'^ and use (p8|). It writes 

[E ^"'"* (E ^^^- c°«(^^^) + ^tfe «i^(^^*)) ) 

X ( E ^"'"* ( E ATfc cos(;0f) + B^, sin(/0t) j j . 



16 RAPHAEL COTE AND STEFAN LE COZ 

Now, the multinomial development gives 

n 



X 

fc=l 



^ (a+j^(x) cos(/6't) + B+i^{x) sm{Wt) 

k 






jfc 



\i=0 



Fix some {ik)k, {jk)k and define the decay rate k = J2k=i ^C**; + J*;)- Then 

'^ ^ JZ^**^ + ■^'^) = n + (m - n) = TO > 2. 
fc=i 

The product factor is a trigonometric polynomial in t, it can be linearized into a sum of sin and 
cos with frequency £9 and i ^ ^j. fc(zfe + j^) — k. 

Of course, as W^ G Jif{R^), the higher order terms (i.e. with k ^ iVo+1) all fit into 0(e"(^°+i)'"). 

It is now clear that Aj,^ and Bj,^ are polynomial in Aj^., Bjj^, Pn.rm and Qn,m- It remains to 
see that the Aj_k or Bjk that intervene (i.e u- +jfe > 0) come with fc ^ k— 1. Let a be the maximal 
index such that ia + ja > 0. Recall ii + ■ ■ ■ + iN„ + ji + ■ ■ ■ + JNa = '^i ^ 2. If ia + ja ^ 2, we have 
2a ^ a{ia + ja) ^ k so that (as k ^ to, ^ 2) a ^ k — 1. If ia + Jq = 1, there exist b ^ 1, b ^ a, such 
that ib + jb ^ ^ and 

K = E /c(ife + jfe) ^ a(ia + ja) + b{ib + jb) ^ a + 1. 

Finally the product has the desired properties. D 

Proof of Proposition \2^. By definition of W, we can compute: 

No I k 

fc=l \j=0 

+ {L^2Bj,k - j6'^j,fc - kpBj^k) sm{j9t) . 

From the computations of Claim p3, it suffices to solve for all ^ j ^ A: ^ A^o 

LR2Aj^k + jOBj^k ~ kpAj,k = Aj,k, 



(29) 

i LR2Bj,k — jdAj,k — kpBj^k = Bj^k 

Obviously, one starts to solve for A; = 1, then from this fc = 2 etc. so that at all stages Aj^k and 
Bj^k are well defined (remark that Aj^i = Bj^i ~ 0). 

We initialized the induction process by setting Ai^i = aYi, Bii = aY2, and Aqi ~ So.i = 0. 
Assume that Aj^k and Bjk arc constructed up to fc ^ fcg — 1 and belong to Jf (R^), we now 
construct Aj^koi ^i,fco for all j ^ fco- By Claim E4l all Aj^ko and Bj^^o are constructed for j ^ fcg 
and belong to .J^{M?). 

Consider now the operator Lj^k,, = Lc2 — {kop + ij9) Id, Lj^ko ■ ^(C^) -^ ^(C^). As e = p + i9 
is an eigenvalue of Lc2 with maximal real part, for all fco ^ 2 and all j, kop + ij9 ^ Sp{L) so 

that Lj^ko is invertible. Let X = Ljl^^{Aj^ko + iBj,k„), and define C := .^«(X) — I ^*^„_x J, 
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■^■"'{^^ ) \ „„ +!,„+ n r, r- ^mi\ 



D :== J&.(X) = " ^_ , so that C,D € ,J^{W) and X = C + iD. Then we compute 

Aj,ko + ^Bj,ko ^ Lj^koiC + iD) 

= Lr2C + iLR2D - kopC - ikoD - ij6C + j0D 
= (Lr2C - kopC + jOD) + i{Lj,2D - jdC - kopD). 

Hence Aj^ko = C and -Bj.^,, = D are solutions to the system (p9[). D 

We now switch back notation from vector valued functions to complex valued functions and 
summarize what we have obtained. We use again the subscript 1. Hence we can consider V\ ", 
U^° defined by 

V,^"{t,x) := e'^'W^°it,x), U^"it,x) := Ri{t,x) + V,^<>{t,x). 

Then we define 

Err^° {t, x) := idtU^" + AU^" + /(C/f °) 

= idtV,''" + AFi^« + /(i?i(t) + V,""') - fiRiit)) 

= zidtVl"" + ^cVl"" - .^ciV,''")) 

= ie''^\dtW^'^'' + LcW^'' - ^ciW^"))- 

By Proposition H, Err^°{t,x) = 0{e-^'^°+^'^P^). Also, from (H) we deduce 

Vi^°{t) = ae^"*r(t) + 0(e~2p*)^ so that for aU s ^ 0, there exists C{No, s) such that 

(30) Vt^O, |K«(i)||H-,(R.)SCa(JVo,s)e-''*. 

3.2. Proofs of Theorems § and |. 

Proof of Theorem 0. Let Nq to be determined later, we do a fixed point around Ui °{t). Suppose 
u = U^°{t) + w{t) (with w{t) -^ as t -^ +oo) is a solution to ( |NLS|) , then 



idtw + Aw + f{U^° + w) - /(C/f ") - Err^° {t) = 
From this, Duhamel's Formula gives, for t ^ s, 

w{s) = e'^^'-'^w{t) + i Te'^^^-") (/((C/f" + w)iT)) - /([/^''(t)) - Err^" {t)'^ dr, 
so that 

e-'^'w{s) = e-*^*ii;(i) + i / e"*^" (/((C/f'" + ?«)(t)) - fiU^^ir)) - Err^° (t)^ dr. 
Letting s — ?► +oo, as 'w{s) -^ 0, we are looking for a solution to the fixed point equation 

<.+oo 

w{t) = -z / e^^(*-^)(/((t/f" + w){t)) - /(t/f"(r)) - Err^«{T))dT. 



Hence, we define the map 

e''^^*~"H/((i?i + l^i"^" + v){t)) - /((i?i + 1/i^")(r)) - ^rrf°(T))dT. 

Fix cr > |, so that i?°'(R'^) is an algebra, and let B,To to be determined later. For 
w e C((To, +00), iJ'^(M'')) define 

||u;||x^_ =supe(^°+i)^*|iz«(i)||^.(M^), 

to be the norm of the Banach space 

]w\\x- ^,_ < +00I. 



Consider the ball of radius B of X^ ^ 



[weC{{To,+oo),H''{R'')) 



18 



RAPHAEL COTE AND STEFAN LE COZ 



By (30), we can assume Tq is large enough so that 

l|Vf'°||if^(Rd) ^ 1 and also Be-'^^°+^^P^'> < 1. 

Our problem is to find a fixed point for ^, we will find it in X^ ^ (B) for adequate parameters. 

Notice that for t ^ Tq, \\V-^^''(t)\\H^(TSid) ^ 1. Hence, we will always work in the H'^ {W^yhall of 
radius r^ — ||'&i||ff"(R<i) + 2. Due to C^^^ smoothness of /, there exists a constant Kg- such that 

Va,6e BH''(Rd){r^), \\f{a) - f{b)\\H-(R-i) < ^a||a - 6||//.(Rd). 

In particular, for all t, 

\\f{Rl{t) + V.^'^'it) +V)- /(i?,i(t) + yi^«(i))||ff.(K.) < K^vllH^^R^y 

Hence, as e*^^*^*) is an isometry in i/'^(M'^), for any v £ X^^ ^ {B) we have 



l*H(i)llH<'(R'i) 



,*A(t~T) 



t 

+ 00 



JVn 



No 



f{Ri + V^'° +v)- f{Ri + V^'°) - Err'^° {T)dT 



H'iR'i) 



< / (||/(i?i + y/^«+^;)-/(i?i + y/^'')||^.(R.) + |l£;rrf"(r)||^.(„.))d- 



No\ 



.No, 



< 



< 



+ 00 



iK,\\v\\H.(R.) + CiNo,a)e-^''"+'^P^)d7 



K,B + C{No,a) 
(No + l)p 



-iNo + l)pt 



First choose Nq large enough so that , . ,"^i . ^ i Then choose B = 2 , ., "fJ . Finally choose Tq 
large enough so that C{No, cr)e~^^° ^ 1. Hence we get 

This shows that ^ maps X^ jy (B) to itself. Let us now show that \I' is a contraction in X^ ^ (B). 
Let v,w G X^ jY (B) then we have 

*(v)(t) - *H(t) = -^ / e"^^'~'^fiRi + V^^° + w) - /(i?i + Vi^° + w))ds. 



As previously, we have 



^{No + l}pt 



+ CX3 



e*^(*-^)(/(i?i + 1^1^° +v)- f{Ri + Fi^° + w))ds 
< e^^o+i)^* / ||/(i?i + Vl"" +v)~ f{R, + Fi^° + u;)|U.(R.)ds 



+ 00 



H^ 



+ 00 



^ ^{No + l)pt / K,\\y(^s) - w{s)\\H.(^.)ds 



+CX3 



Jt 



II ll-^To.JVo 



-{No + l)pt 



€ 



K„ 



, ■^\\v — w\\X'' 

(A^o + l)p "■«•"» 

Taking the supremum over t ^ Tqi we deduce that 



||*(w) -*H|| 



X, 



^- ■:—\\V^W\\X'' ^ —\\V ~ WVxX" 
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Hence, \[' is a contraction on X^ ^ (i3), and has a unique fixed point v. Notice that we have 
obtained a unique fixed point for any a '^ ^: from this we deduce that v does not depend on a, 
and hence, v G C°°{[To,+oo),H°°{R'^)). Then u = Ri + V^^" + 5 is the desired sohition. D 

Proof of Theorem H. The proof is essentially a generalization of that of Theorem 0. Let ut, to be 
fixed later and assume that Vi, > v\f. Let A^o to determined later and a e M, from this we dispose of 
a profile V^ °{t), Ui "(t), an error term Err^ °{t) associated to Ri{t), and an eigenvalue X = p + i6 
of L. We look for a solution of the form u{t) — U^ °{t) + X]j>2 ^j(*) + wit). Then w satisfies 



idtw + Aw + f{U^° + J2 R] 

J>2 



w) 



f{U^»)-J2f{R,)-Err^« 

J>2 



0. 



Hence considering the map 

r'+oo 



V !->• *(«) 



J>2 j>2 



we are looking for a fixed point for ^, in the set X^ ^ (B) (defined in the proof of Theorem g) 
for adequate parameters TQ,NQ,B,a. Let cr > |. As previously, let Tq large enough so that 
\\V^°{t)\\H-{M-i) ^ 1 for t ^ To, and Se"(^«+i)''^o ^ 1, so that we remain in a baU of radius 1 in 

Using exponential localization of the solitons Rj and of the profile Ui " , we deduce as in the 
proof of Theorem ||that for some Ka = K{f, \\Ui °||//^(Rd) + Ei>2 \\Rj\\H'{m^) + 1), we have 



l/((c/f« + J2Rj+ «)) - f(u^°) - E fiRMH^m ^ K,\\v\\„.^t^^) + o{e 

i>2 ]^2 



— 2a^/uj^Vi:t 



), 



possibly by taking a smaller value of uji, such that cj^ ^ ai, where ai is the (exponential) decay 
rate of Ui ° . Notice that ai is independent of A^o, due to the construction of Ui " . Hence we have 
as in Theorem 0: 



|*(t;)(t)||^.(Rd) 

f + 00 



/ 



JA{t-s) 



< 



€ 



t 

+ 00 



(/(C/f ° + ^ i?,. + „) - /(C/f °) - ^ f(R^) - Err^'>)ds, 



J>2 



j>2 



H" 



\I{ui 



No 



J2R,+v)-f{U^")-J2f{R,)\\H' 

j>2 j^2 



\Err^nH^ 



4-00 



(if,||t;||ff.(R.) + C(iVo,a)e-(^«+i)''^ + C{a)i 



— 2a.-^/u)^v^s 



)ds 



K^BjrC{No2^ (No + l)pt , ^(^) „-2a,/ZJ7f.t 



(A^o + l)p 



2aJuuv-i, 



First choose Aq large enough so that ,f/^Zi\ ^ \ £^nd set S := yiv+T) ■ Recall that w* > wt|. We 
chose v\^ large enough so that from the choice of uj-^^v^,, we have 

-=^ < — , and 2a./uJlVi, ^ (Nq + l)p. 

2aJu^Vi, 3 



Finally choose Tq large enough so that 

Be-^'^o+^^P^o < 1, and C{No,<T)e-P^° < 1. 

From this, ||*(z;)(i)||ff<,(Rd) s^ Be~'^^°+'^^P* for t ^ Tq, i.e. ^ maps X^^j^^{B) to itself. 
Similar computations show that ^ is a contracting map, so that it has a unique fixed point 
w. Again as in Theorem g, u) does not depend on a and w G C°°([To, +(X)),iJ°°(M'')). Then 
u = C/-^ ° + X]i>2 Rji^) + ™(^) fulfills the requirements. D 
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Appendix A. Exponential decay of eigenfunctions to matrix Schrodinger 

OPERATORS 

We consider an operator L : H^{R'^,C^) C L^{R'^X^) ^ L^(R'^,C^) of the form 

Wi -A + w + Vi" 



^ \A-uj + V2 W2 

where a; > and Vi, V2, Wi, W2 are complex- valued potentials satisfying the following assumptions. 

(VWl) There exists q e (max{2, f }, +00] such that Vk, Wk e L«(R'*) for fc = 1, 2. 
(VW2) lim|,|^+o, Vk{x) = lim|,|_^+oo Wk{x) = for fc = 1, 2. 

Assumptions (VW1)-(VW2) are probably not optimal, but they are sufficient in the context in 
which we want to apply the following Proposition Eq. 

Our goal is to prove that if L has an eigenvalue which does not belong to the set 
{iy, y e R, |y| ^ 1^} (which is the essential spectrum of L, see e.g. Q) then the correspond- 
ing eigenvectors are exponentially decaying at infinity. Note that it was previously known only 
for eigenfunctions corresponding to eigenvalues lying in the strip {z e C, |jS»(z)| < w} and with a 
restricted class of potentials (see [^ ) . 

Proposition 25. Assume that (VW1)-(VW2) hold. 



Take u,v Cz H {W ,<C), A G C \ {iy,y £ M, \y\ ^ uj}, and suppose that for U :— { \ we have 

LU — XU . Then there exist C > and a > such that for all x G M'' we have 

\u{x)\ + \v{x)\ s$Ce-"l^l. 

Our proof consists in obtaining estimates on fundamental solutions to Hclmholtz equations and 
considering the eigenvalue problem LU = XU as an inhomogeneous problem. 

A.l. Fundamental solutions. For a given /i e C, a fundamental solution of the Helmholtz 
equation in M'^ is a solution of 

(-A - Ai).g^ - ^0. 
Setting I' := ^^ fundamental solutions of the Helmholtz equation are given by 

2\x\'^[2tt)2 

where Hi is the first Hankel function (see e.g. [Q). For fj, = pe*^ with p ^ and 9 G [0,27r) 
we defined y/JI by y/JI :— p^e^^. Defining ^/^ in this way ensures in particular that g^ is square 
integrable for /i ^ M+ . The fundamental solutions gf^ verify the recurrence relation 



9 „d 



^ ^ '~ 2n\x\ ■ 



We deduce the following formula for the fundamental solution. For d = j + 21 where j = 1,2 and 
I €N\ {0}, we have 



(31) 5^+''=E«'(-1)'(5m) 



(fe)l^|-2;+fc 



k=l 



where the coefficients (af) are positive and the exponent (k) denotes the fc*^ derivative. 

Lemma 26 (Estimates on fundamental solutions). Let /i G C \ IR+. Then there exists t > and 
C > such that 

\gtix)\^Cgt,ix)forallxeR''\{0}. 

In particular, gf^ is exponentially decaying at infinity with decay rate ^/t, i.e. |gf (2:)] ^ Ce~^'^' 
for \x\ large enough. 

We separated the proof of Lemma Eq into two proofs depending on the oddness of d. 
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Proof for odd d. We have y^ — pae*^ . Choose r > such that ^/t = p^ sm |. It is well-known 
that gjtix) = Y7=e^^^^y It follows from easy computations that 

Since 



2ypsin| 

this readily implies that for all a; G K'' we have 

\gl,ix)\ ^ CgUx), 

which proves the lemma for d = 1. 
Similar calculations lead to 

(32) KgMWKCC-lflsijC^^forallfce 



Assume now that d ^ 3 and take Z G N \ {0} such that d = 1 + 2L Combining (p\\) and (32) gives 



\gl+''{x)\ ^ Cg]X^\x) for all x e M"^ \ {0}, 
which is the desired conclusion. D 

Proof for even d. Let i/ G N and z G C We have the following asymptotic expansions on the 
Hankel functions (see m). 



iHliz) R 


i--ln(z) 

TT 




for \z\ close to 0, 


iHl{z) « 


v\z-'' 




for \z\ close to 0, i^ 7^ 0, 


2-''TT 


Hliz)^ 


M/^e*- 


^-f) 


for \z\ close to + 00. 



Therefore, we can infer the following estimates on the fundamental solutions. Recall that d — 2 + 2v 
and p = pe'^. 

(33) \gl{x)\ ^ C\ ln(p^|a;|)| for |a;| close to 0, 

(34) \g'l{x)\ ^ C\x\-'' for |a;| close to 0,z/ 7^0, 

(35) \g'l{x)\ < C\x\''^''+^^e'P^ ^^'^(i)M for |x| close to + 00. 
For T > 0, the function gl^ verifies gl^ > and 

(36) g-Ax) w C\ ln(r5 |a;|)| for \x\ close to 0, 

(37) gtAx) ~ C\x\-'' for |a;| close to 0, i^ 7^ 

(38) gt^x) « aixp^'^+i^e-^'l^l for \x\ close to + 00. 



Choose T > such that t^ = .ypsin|. Then we infer from (|35|)-(p8D and the continuity of 
fundamental solutions that there exists (7 > such that 

\gi{x)\^CgiAx)ioT^nxeR^\{0}, 

which finishes the proof. D 

A. 2. Exponential decay. We start with a regularity result on eigenfunctions. 
Lemma 27. Assume that (VWl) is satisfied. Take A G C \ {iy, y G R, \y\ ^ lo}, u,v e H'^(R'^, C) 
and assume that for U :— I \ we have LU = XU . Then u,v G W^'^'^(R'') for any r G [2,q]. In 
particular, u,v E C^{M.'^) and \iin\x\_^^rx> u{x) = lim|^|_j.^oc v{x) = 0. 
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Proof. The result follows from a classical bootstrap argument. Let the sequence (r„) be defined 

by 

f ro = 2, 

< 1 _ 1 I d-2rj 



{ rj + i q dvj ' 

where q is given by (VWl). An elementary analysis of (r^) shows that there exists jo such that for 
all < j < jo we have Tj+i > rj, -^j^ > and ~^J'° < 0. 

By induction, it is easy to see that for all j — 0, ..., jo we have u,v £ W^'^^(R'^). For j = it is 
by definition of u, v. Take any ^ j < jo and assume that u,v G W'^'^^ (R'^) . Since ~^''^ > 0, by 

Sobolev embeddings we infer that u,v £ L'' ^''j (M ). Then, (VWl) and Holder inequality imply 

Wiu, Viv, V2U, Wiv G ^'■^■+1 (M'*). 
Combined with U = (u, v^ satisfying LU = XU, this leads tou,v e Ty^.rj+i j^j^d^^ 

In particular, we have u.v G W'^'^^oiM.'^). Since -—. — — < 0, from Sobolev embeddings we infer 
u,vG L°°{R'^). Then by (VWl) we get 

Wiu, Viv, V2U, Wiv e L^R''). 

As before, combined with LU = \U, this leads to u,v G W'^''^{R'^). The conclusion follows by 
interpolation. D 

For the rest of the proof, it is easier to work with the operator 

y W2 A-uj + V2 

where P = (\ ^ ■ ) ■ The potentials V{,V^,W{,W^ verify also (VW1)-VW2). The spectrum of 
L' is Sp(L') = Sp{iPLP^^) = iSp(L). Hence if A e C is an eigenvalue of L with eigenvector U 
then A' := iX is an eigenvalue of L' with eigenvector U' — I / I ■= PU. 
Write L' - X'l = H + K where 



'2 ^2 



Take 






It is well known that we can represent u' and v' in the following way 

u' = gt^+x' * h and v' = -gt^^x' * /2- 

Let /ii := — w + A' and /i2 := ^'^ — X' . From the assumptions on A' we infer that //i,/i2 satisfy 
the hypothesis of Lemma EO. Let ri, T2 be given by Lemma Eq and set r := minJTi, T2}. Take 

^ ._ rh\ . (\h 



J2J ■ V1/21 

u := g'Lr * h and v := gt^- * f2- 
Claim 28. There exists C > such that 

\u'\ ^ Cu and \v'\ ^ Cv. 
Proof. This readily follows from Lemma E6l D 

Lemma 29. Set w := u + v . There exists C > and a > such that 

w{x) s: Ce""l"^l for all x e R"*. 
The proof of Lemma E9 follows closely the proof of Theorem 1.1 in Km. 
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Proof. Set / := /i + f-z- We first note that w E C'^(M'*). Indeed, by definition w satisfies 

(39) - Aw + rw = /. 

Since, by (VWl) and Lemma |^, / G L«(M'*) this implies u; e 1^^'''(M'^) and in particular 
w G C°(R''). 

Now, we claim that there exists R> such that for all a; G M'* verifying |a;| > R we have 

(40) rw{x)-m r 
^ ' w{x) 2 

Indeed, setting T{x) := {\V{\ + \V^\ + \W{\ + \W^\), by Claim || we have 

/ s$ T{x){\u'\ + \v'\) s^ CT{x){u + i) = CT{x)w. 

Therefore 

TWJX) - f{x) 

-— ^ T - CT{x). 

w[x) 

By (VW2), we can take R large enough so that CT{x) ^ -^ for |a;| > _R, which proves (pol). 

Note that w > by definition. Since w G C"(R'^) n VK2,9(]gd)^ i^^ere exists C^ such that for all 

X eW^ with |x| < Rwe have 

SC u;(a;) sC Cr. 

Define ^(a;) := C^e^^^l^l'-^^ It is easy to see that 

(41) -AV' + ^V ^ on R'^ \ {0}, 

w;(a:) - ip{x) ^ on {x G R'^, |a;| < R}. 

Therefore we only have to prove that w{x) ^ ip{x) for |x| > R. We proceed by contradiction. 
Assume that there exists xq G R'^ with |a;o| > R such that w{xo) > V'(a^o)- Define the set 

n :={a;GR'',-u;(a;) > V^a^)}- 

Then 51 is a non-empty open set, for all x G Jl we have |x| > R and for all x G dfl we have 
w{x) - ipix) = 0. On n, by (H), ^ and (|l|) we have 



A(u) — ^) = Aw — A?/) — Tw — / — Ai/) 

TW — f T 

= -w - Ai/j > -(w - -0) > 0. 

w 2 

By the maximum principle, this implies that w — ip ^ on il, a. contradiction. Thus, for all x G M*^ 
we have 

w{x) s; tP{x) = Cfle^VId^l-^) ^ ^^eyT«e-y¥kl ^ Ce~VTl^l. 

This ends the proof. D 

Proof of Proposition ^5. The statement is an immediate consequence of Lemma UTl Claim ES and 
Lemma E9^. D 

A. 3. Higher regularity and decay. Upon assuming more regularity and decay, we can obtain 
more regularity and decay on the solutions to (L — XI) = A. 
The new assumption is the following. 

(VW3) Vi,V2, Wi,W2e .Jf{C). 

Recall that ^ was defined in (EST). 



Proposition 30. Assume that (VW1)-(VW3) hold. 

(i) Let A, u and v he as in Proposition W5. Then u,v € J^{C). 

(ii) Let X ^ Sp(L) and take A G ^(C^). Then there exists X G ^(C^) such that 
{L- Xld)X = A. 
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Proof, (i) The assertion follows from similar arguments to those used in the proof of Proposition [ 



provided we remark that (using the same notations) D'^u' = g'^^<x'*D"^fi, D°-v' = —9'^i^^x'*D"-f2 



and D°-fi,D°-f2 satisfy the same properties as /i and /2. 

(ii) Since A ^ Sp(L) the operator L — Aid is invertible, hence the existence of X g if^(M'^,C^) 
such that [L — AId)X = A. Regularity of X follows from a standard bootstrap argument as 
explained in the proof of Proposition p^ (ii). We now recall that L = —iP^^L'P. Hence, if we 
define X' = PX, A' = iX, and A' = iPA then 

{L' ~ X' ld)X' ^ A' . 
Recall that L' - X'l = H + K. Set F = ( ^M := KX' and A' = ( "M . Then we can represent 
X' = { ] in the following way 

Xi = 9-^+x' * (yi + ai) and X2 = -Q-^-y * (^2 + 02)- 

The terms g'^^.y * ^1 and g''i^_y * ^2 are clearly exponentially decaying, with decay rate a. Since 
Vi, V2, Wi, W2 S =^(C^), it follows that each component of Y is also exponentially decaying with 
rate a. Hence g'L^,x' *2^i and g'L^^y *y2 are exponentially decaying with decay rate a. The decay 
rate of the derivatives of X' is follows immediately if we remark that for any multiindex a we have 
D^Xk^gt^^y*D-{yk+ak)iovk = l,2. ' D 

Appendix B. Instability of solitons and multi-solitons 



Since ( NLS ) is Galilean invariant, we can assume in this section without loss of generality that 
ui = a;i = 71 = 0. Hence Ri{t,x) = e*'^i*$i(x). 



Recall that, as defined in Section 3.1, Y{t) is of the form e '^*{cos{0t)Yi{x) + sin{9t)Y2{x)), 



where Yi,Y2 are smooth, exponentially decaying functions, along with their derivatives. No- 



tice that if u{t,x) is a solution to ( [NLSD and T,'d G R, then so is ^(r - t,x)e'^. The 
hypotheses of Theorem g are verified by $1 and therefore also by $1. Hence the conclu- 
sion of Theorem | holds for Ri{t,x) := i?i(-i,x) = e'"i*$i. Let u G 'Tf{[To,oo),H'^(R'^) 
be the solution constructed in Theorem g associated with the soliton Ri{t,x) and correction 
e"''*(cos(6'i)Yi(a;) -I- sm{et)Y2{x)) + 0{e~^P*) (i.e. u = wi in the notations of Theorem ||). In 
particular, for all cr 5^ 0, 

Vt ^ To, llu(t) - Ri{t) - r(t)|la.(R<t) «; Ce-^p\ 
Note that we construct u on Ri and not i?i so as to have instability forward in time. 

B.l. Orbital instability of one soliton. First let us prove a modulation lemma. 
Lemma 31. There exist e > 0, Iq ^ Tq and M ^ such that 

inf ||u(to) - $i(x - y)e*''||i2(B(o,M)) = e > 0. 

Proof. Let io > ^0 to be determined later. Up to increasing tg, we can assume that wiip = 0{2tt). 

Consider Q{y, •&) = ||u(io) — $i(x — 2/)e"'||i2(]jd). The function 9 is continuous on IR'^+^. Notice 
that for ?9 = and y^O, one gets 9(0, 0) < Ce"''*" . 

Now, we have that liminf|j,|_j.oo inf,jgB9(j/,i?) ^ 2||l'i||^2(]{d) — Ce^''*" due to space localization 
of $1, so that, as "i? e R/27rZ compact, if to is large enough, infj,gRd ,jgj{ Q(jJ, "&) is attained at some 
point (yo,?9o). 

Assume 9(yo,i9o) = 0, i.e. u(io) = ^i{x - yo)e'^°. 
Claim: There exists a continuous function 77 such that 77(0) = and |yo| + |7?o| ^ r]{e~P^°). 

Indeed, first consider y^. Denote g{y) = |||$i| — |$i(- — y)|||^2™d)- We have 

= 9(yo,^o) ^ |||u(io)| - |^l(- - yo)\\\L-(R^) > lll^ll - l*l(- - 2/o)|||L2(Rd) - q|y(to)||L^(R.). 

As ||y(io)||L2(ud) < Ce^P*", we get 5(2/0) ^ C^e^^^*°. Now, due to space localization of $1, 
g{y) -> 2||$i||22^jj^^ > as \y\ -^ +00. Let (y„) be such that g(y„) -^ 0, and y„ 7^ 0. Then up to 
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a subsequence, yn — >■ y°° and g{y°°) — 0, so that |<i>i| is periodic and as $i e L^(R'^), $i = 0, a 
contradiction. This shows that y — ?► as ^(j/) — > 0, and it gives the bound on j/o- For i?o, 

= ||u(to)-$i(--yo)||L2(Rd) ^ -||u(to)-$i||L2(R<i) + ||l>i-l>ie'''«||i2(Kd)-||l>i-$i(--yo)||L2(R<i), 

As ||$i - $ie''''«||i2(R.) = |1 - e*''«|||$i|U2(R.), we deduce that \M ^ Ce-P'- + Cg{yo). This 
concludes the proof of the claim. 

Denote Tj,^^ the tangent space of ^ := {l>i(- - y)e'^\{y,-d) £ R'^} at point $i. Note 
that, due to the Claim, =^ is a manifold. It is easy to see that T^^^ C keric (by differ- 
entiating the relation A$i(a; — y) + g{\^i{x — y)\^)^i{x — y) = wi$i(a; — y)). But for all t, 
{cos{9t)Yi{x) + shi{9t)Y2{x)) ^ keric (as 11,1^2 are build on an eigenvector for an eigenvalue of 
positive real part of Lc). As u(io) = $i + eP*-» {cos{eto)Yi{x) + sin{eta)Y2{x)) + ©(e'^''*"), up to 
choosing to + 2fc7r/0, (fc G N large) instead of io, this proves that u(to) ^ ^- We proved that for 
io large enough, 

inf ||u(io) - ^i{x ~ y)e'^\\m{9.<i) > 0. 



Assume that this does not hold when we restrict to L^{B{0,M)), for any large M. This would 
mean that for all to ^ 0, there exist y^ S K'', ^m G R such that 



||u(to) - *i(a;-y,„)e*''""-||i2(B(o^„j)) 



€ 



Then by localization arguments, (?/m) remains bounded, so that up to a subsequence, y^. — > J/oo 
t^m -> t^oo- Therefore ||u(io,a;) - $i(x - 2/oo)e"'°°|lL2(Rd) = 0, so that u{to,x) = $i(a; - yoo)e*''°°, 
a contradiction. D 

Proof of Corollary IJ. Let tg and e be given by Lemma ^ Take an increasing sequence (S'„) so 
that S'n ^^ +00 as n — >■ +00, and define r„ :— S*,! — ig and 

Unit,x) ■.^uiSn-t,x)e^"^'^-. 

Then w„ e '^([0,T„],iJi(K'')) is a solution of ( |NLSD , and 



M„(0, x) = u(5„, x)e-*'^i^" = $i(a;) + O^. (e""^"), 
Un{Tn,x)^u{to,x)e~'"''^-. 

Therefore, ||wn(0) — Ri{0)\\H'(m'i) ^ Ce^P^" -> as n — ?> +00. Due to Lemma ^, we deduce that 
for all n S N we have 

inf \\uniT„) - e'''$i(- - y)\\ms.^) > inf ||u(to) - $i(x - y)e'^\\LHB(o,M)) ^ e, 
yeR'',i?eR ^ ' yeW^^em k \ • n 

which is the desired conclusion. D 

B.2. Instability of multi-solitons. 

Proof of Corollary 0. Let T > 0, M be given by Lemma ^ and e, (u„), (r„) be given by Corollary 

I 

The idea is the following. We use the fact that u„(r„) is e-away from the orbit of the soliton 

Ri. Given a parameter /, we consider at time / an initial data w{I) which is Wn(0) adequately 
shifted, denoted by u„(/), plus the sum of the Rj{I), j ^ 2. (All the functions will depend on n 
and /, although we do not always show this dependence for convenience in the notation). We aim 
at controlling w up to time / + r„. The role of / is to ensure that the interaction of u„ and the Rj 
are small: as {un{t)\t e [0,r„]} is compact and the Rj{t) {j ^ 2) are localized away from u„(t), 
their interaction goes to as / — > +00. Using a Gronwall type argument, we are able to show that 
w{I + T„) is Un{I + Tn) + Z]j=2 -^j (^ + ^") + o/^+oo(l)- As u„(T„) is £-away from the soliton 
family, we deduce that w^I + Tn) is e — o/_^-|_oo(l) ^ £/2 away from the family of a sum of solitons. 
Given I ^T, define u G "^([7,/ + TJ, iJi(M'')) by 

Un{t,x) = Un{t - I,x). 

Possibly increasing / so that ujil = 0(27r), 

we have ||u„(/) - Ri{I)\\H<'{Rd) = ||wn(0) - RiiO)\\H'{R'i) -^- as n ^ +00 and u„(/ + r„) is 



26 RAPHAEL COTE AND STEFAN LE COZ 



e-away from the $i-soliton family. Consider the solution w„ S '^([/,r*),i?^(M'*)) to ( |NLS[ ) with 
initial data at time / 

N 
i=2 



If T* < +00, the blow-up alternative for (NLS) automatically implies instability on the 



multi-soliton, hence we assume T* = -1-cx). Let a > d/2 be an integer. Notice that, as 
u„ e '^([0,r„],iJ'"(R'')) and [0,T„] is compact, the set {M„(i)|i G [0,r„]} is compact in H''(R'^). 
In particular, supjgroj-^i \\un{t)\\H''{\x\^R) — J' as i? — > +00. Hence, as the Rj are decoupling as 
time grows, there exists a function 77(7) such that 77(1) — > as / — > +00 and 

Vte [/,/ + r„], ^||u„(i)i?j(i)||ff" <ry(/). 

J>2 



Denote Xj(t) = Vjt + Xj. Up to modifying the function 77, we can also assume that the Rj{t), j ^ 2, 
are far away from Xi{t) = 0, and that the multisoliton R{t) is near the sum of solitons X]i=i ^ji^)-: 
that is 

AT 

Vt ^ /, J2 \\Rjit)\\HHB{0.2M)) + \\R{t) - ^i?j(i)||H-(K.) S; Vil)- 
j=2 i=l 

Finally we denote J = I + Tn and 

AT 
Z{t) = Wn{t) - {Un(t) + ^i?j(t)). 
J=2 

Now, as / is '^°°, for all i? > 0, there exists C{R) such that 

(42) Va,6ei3(0,i?), \f{a + h)-f{a)-f{b)\^C{R)\a\\b\. 

Indeed, this expression is symmetric in a, h, so that we can assume without loss of generality that 
l&l s; \a\. As /(O) = /'(O) = 0, we have that |/(&)| < C\b\^ ^ C\a\\h\, and a Taylor expansion shows 
that 

\f{a + b)-f{a)\^b f \f'{a + tb)\dt^b sup \f{x)\^C\b\{\b\ + \a\)^C\a\\b\. 

Jo a;GB(0,|a|+|f)|) 



Now, as H'^{M.'^) is an algebra, we deduce from (42) that there exists a constant C > (depending 
only on the $j) such that for t € [/, J], 

N N N 

\\f{z{t)+Un{t)+Y,Rj{t))-f{Ur.)-Y.f{Rjm\H^m^C\\z{t)\\H^^m^)+CY,\\^n{t)RjmH^{m^)- 
j=2 j=2 j=2 

The function z satisfies the equation 

ZZt + AZ + / ( 2 + 7i„ + ^ i?, j - /(U„) - ;^ f{Rj) = 0, 

Since z(/) = 0, Duhamel formula for z gives 

z(t) = II e^^(*-^) (f (^z{s) + u,M + f2 R^i'-^)) - /(""G'^)) - E fiRji'))] ds- 

Hence, for all t E [I, J] 

IkWIIff-HR'') ^ C* / (lk(5)llH-(R'i) + '7)c^s < C" / ||z(s)||^.(Rd)ds + ?7(/)(t - /). 
By Gronwall's Lemma, we deduce that for t £ [I, J], we have 

||z(t)|U.™.) < C77(/)(t - /)e^(*-^) < C'„77(/), 
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where C„ = CTnC^'^'^ . Thus for aU n e N we have 



N 



Un{J) ~ Un{Tn) ~ ^ -Rj(>^) 
J=2 



^ CnTlil)- 



//"(R-i) 



Now choose /„ such that C„77(/„) < e/3 and set J„ = /„ + T„. Then ||z(J„)||^<T(]{d) < e/3. Then, 
given yj e K'^, i?j G M, we have (denote Cj = Cj(i) = —^bjP* + ^ji + 7o) 

TV 



3 = i 



> 



> 



N 



N 



ln{Tn) + J2 R'J^Jr.) - J] $,(• " yM^"^^^"" 



+'»j) 



3=2 



i=i 



L2 



TV 



WniJn) - Mn(T„) - ^ RjiJn] 



i=2 



L2 



AT 



,(r„) ~^,{x- yOe'-'i + ^($,(x - x,(J„))e^(^^-+^^) - $,(x - y,)e*(^^-+''^)) 

i=2 



r/3. 



L2 



Now consider yj, dj that reahze a near infimum, say ||w„(J„) — ^ ■ -^ $j(- — yj)e*'^2^j'^+''j)||^2 ^ 2e. 
Then considering the L^ norm on baUs B{xj{Jn), R) around each exited state Rj, j ^ 2 (for some 
large and fixed radius R), we see that, up to a permutation if two $j or more are equal, we must 
have yj — Xj(Jn) — 0(1) for j ^ 2. In particular, this implies that 



N 

^($j(a; - a;j(J„))e'(5''^-^+^^) - ^j{x - ?;j)e'(^^-^+''^) 

J=2 

up to increasing again /„. Thus, 

N 



H'{B{0,M)) 



N 






o/„^+oo(l) sS e/3. 



iVj-X + ^-i)! 



j = l,...,Ar 



J=l 



i=i 



AT 



2^i?j(/„) - i?(/„)||/i-<T(Rd) -> 0, 



> |lu„(r„) - $i(a; - yi)e^''i |1l^(b(o,m)) - 2£/3 
^ e - 2£/3 ^ e/3, 
where we used Corollary on the last line. As 

||'(«n(-?^n) - R{In)\\H''{R'^) =^ ||w„(/„) - ^^ i?j (/„) || /j^ (Rd- 

w„, /„ and Jn satisfies the conditions of Corollary El D 

Remark 32. Notice that we did not use any high speed condition on the Vj . The most delicate point 
here is that we have no uniform spatial decay on Un (as well as on the multi-soliton constructed in 
Theorem W , apart that coming from H°'{M.'^) compactness. We conjecture it should be exponentially 
decaying (in space) around the soliton (resp. every soliton Rj); a proof of this should be related to 
uniqueness of the multi-soliton in the L^ sub-critical case, which is currently an open problem. 

Appendix C. Coercivity for a soliton 
This section is devoted to the proof of Lemma [I4 
Proof of Lemma \lSi. We first remark that Rq is solution of 



(43) 



Ai?o + ojo 



M 



i?o-/(i?o)+i«oVi?o-0. 



Therefore it is a critical point of the functional Sq defined for w G H (Mr) by 



So{w) := :r||Vw||i2(Rd) + - Wo 



l^'ol 



lL2( 



F(w)dx — -Wo • -^'-^ / wWwdx. 
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The quadratic form Hq is precisely 

Ho{t,w)^(^S'^{Ro)w,wy 
Consider z such that w ~ q-^(2^o-x-:^\vo\ t+uiot+io) ^^^^ _^ ^^^ _|_ ^^y Then it is easy to see that 



It is well-known that up to a finite number of non-positive directions Ho{z) controls the H^(R'^)- 
norm of z. Indeed, the self- adjoint operator corresponding to the quadratic form Hq (viewed 

on i/^(]R'',]R^)) is a compact perturbation of I . I, hence its spectrum lies on 



-A + wo^ 

the real line and its essential spectrum is [luq, +00). Since in addition the quadratic form Hq is 
bounded from below on the unit _L^(]R'')-sphcrc, the corresponding operator admits only a finite 
number of eigenvalues in (— oo,a;g) for any cjq < ujq. In particular, there exist Kq > 0, z^o G N and 
X^,...,X^° G L2(]grf) gych that ||X*=||i2(Rd) = 1 for any k and 



2 



3 

2 


l^^llL2(E<i) + 


3 


M\, 1,2 

4 W^'WlH 


T.J 


•^0 

4- .„^ 1 T^^ \~ 


' / 


„„ Vkli-W"^ 



k=l 

Since 

there exists Kq > such that 

lkll^i(R.) ^ KoHoit, w)+Ko}2 (^' ^0 W)l2(R.) , 

where X^{t) := e*(3"o-=^-3l«ol'*+"ot+7o)xfe(x - vot - xq). □ 
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